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NEUTRON SLOWING DOWN IN 
HYDROGENEOUS UEDIA 

L*N« Yurova, A.A* PolJakoT, S.B* Stepanov, G.A. Alimov 

Fast neutron slowing down and the process of thermal ization 
In hydrogeneous media wez*e thoroughly studied in many papers. 

But still it is difficult to say that problem has been studied 
to the bottom. For instance, numerous experimental results of 
investigations of the neutron slowing down In water-uranium 
media do not agree rrith the calculated results. There is limited 
information of the Influence of the medium heterogeneity on the 
neutron slowing down, the thermalizatlon is not studied enough 
in various hydrogeneous media at high temperatures. That is why 
the new ozperimental information is of great selentifio and prac- 
tical interest for understanding the physics of the reactors with 
a hydrogen moderator. 

This paper deals with the results of the experimental deter- 
mination of two Integral values In various media, lengths of 
fast neutron slowing down and diffusion coefficient of thermal 
neutrons . 

I, FAST NEUTRON SLOWING DOWN 
BT-pftT»4m«^ ntal study of the fmst neutron 
slowing down in non-mult inlying media 
^ast neutron slowing down was studied In hydrogeneous media: 
light water (HgO) diphenyl (OjgHjo) isopropyldiphenyl 
as well as In the media ferrum water; alluminium -t- water; 
ferrum + diphenyl and alluminium + diphsnyl. Space distribution 
measurements of slowing down neutrons were carried out in the 
conditions of "limited medium". The term "limited medlm" TnanTin 
the investigated medium of limited size surrounded by an Infini- 
te "reflector". Such geometry of the experiment gave the 
experimental distribution of the slowing dowm neutrons correspon- 
ding to the infinite media with limited quantity of Investiga- 
ted subatances, A similar principle of measurement ?raa applied to 
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tlM itudjing of plijaieal paraaottrs llaltad aecordlog to the size 
of DttltlplylBg aedla placed la the erltloal syatiM. 

A aatheaatlea]. method published la the paper was used for 
a theoretical argument of using the oxperimental method emd for 
the oholee of the system size* / 3 /« As it was expected the 
reaulta of the calculations prored that in the conditions of the 
cylipdrloal geometry of the investigated aedlum placed into the 
infinite "reflector” the diameter of the cylinder must be not 
leas than four lengths of slowing down* (4 ). 

Cylindrical geometry with il cm. In diameter and IIO cm. high 
placed into the infinite graphite "reflector* was used for all 
experimental investigations of the neutron slowing down in non- 
mnltlplylng media* In order to confirm the given geometry 
experimentally space distribution measurements of neutron (Po -f 
Be) source ware made in the following media: diphenyl + graphite 
"reflector" and diphenyl + diphenyl "reflector". The measured 
density distributions of slowing down neutrons up to indium 
resonance In this conditions appear to be identical within 
sxperimental errors. The values of the square length of Blowing 
dowm calculated according to this distributions are correspond- 
ingly equal to: (102,3 “ 3,5) cm^ - for the diphenyl + graphite 
medium said (102,5 - i,l) eo^ - for the diphenyl + diphenyl 
medium. 

Investigating fast - neutron slowing down in aetsQ-hydro- 
geneous media the experiments^, systems consisted of a moderator 
with evenly distributed meted rods. Two groups of media were 
studied: quasihomogeneous media In which metal rods were of 10 mm 
diameter and heterogeneous media with 40 mm rods. 

Measurements were made for neutrons of fission spootnim U-235 
and for neutrons of ( Po + Be ) source, Heutron fisslc n source 
was a convertor made of metallic highly enriched uranium placed 
into a thermallzlng column of a heavy water reactor of the 
Academy of Scienoes of the USSR, 

The apace distributions were measured by 90 mg/cm^ thick in- 
dium foil, the sizes of which were chosen in order to preserve 
the "point" sizes of the source and detectors. Corrections wore 

introduced into the calculated values of the foil activities, 

3 5 6 
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These corrections consemed: !• the oontrlbntlon of the neutrons 
with energies above the first resonanoe level of iiidluB, and 2* 
the corrections concerning the thickness of the neutron fission 
source, ill the measured distributions were calculated to 0,3 mm 
/5/ thickness of the source. Square length of slowing down were 
calculated from the following equation: 



y' 2 y • j ^ "I -i- - 




r 


k*'^^L >x}o(z 




where 


k. X 


- Indium foil activity; 


- constants measured In the experimental neutron 
distribution for an region distant of the source 
where a linear law is present 

The restilts of the square length of slowing down values are 
given In Tables I, II and Figures I, 2. 

The experimental study of the Heutron 
slowing down In a critical system 
In the experimental study of neutron slowing down in a critical 
system was studied In the water-uranium and uranium - Isopropyl 
diphenyl media. Due to neutron leakage measurements, corresponding 
to various geometries of the core. It was possible to determine 
neutron migration areas /6/. The periods of the doubled power sys- 
tem were measured In the experiment with the change of the core 
critical height. According to these periods the values of the 
reactivity were calculated using the formula of inhour. /7/. 

The contribution of delayed neutrons of various energetic groups 
into measured values of reactivity were taken into account. 

values were used in the "inhour*’ formula instead of the 
constant outputs of the delayed neutrons. These values ’^®re 

determined as follows: .zi,x. 

.ef 


where 
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L 


iJe 
L’ Si 


Ti 
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- square length of slowing down of prompt neutron 
fission; 

- square length of slowing down of the delayed neutron 
"1" group; 

- buckling a 


Approved For Release 2009/08/26 : CIA-RDP88-00904R0001 0011 0004-5 





^pprove^o^etease2009^08£6^CIA-RDP88-00904R00010011000^ 



From 




the given values of the reactivity system the 


relations were determined according to the following equation: 


■■ ' 

where H > critical height; 

h - the exceeding of the criti.cal height; 

neutron multiplication coefficient for the 
Infinite media. 

The necessary buckling ralues for calculations were obtained 
from the critical experiments from the measurements of the axial 
and radial density neutron dlstrlbutlona In the core. These dis- 
tributions were measured with copper foils. 

The experimental plant was a zero power reactor, the core of 
which was located In a hexagonal tank male of stainless steel, 
Its side being 50 cm. The core was surrounded by an infinite 
side reflector made of the moderator; there was no upper buttend 
reflector, the bottom of the tank 15 nun being the lower butt- 
end reflector. Fuel elements were hollow cylinders of stainless 
steel filled with (the size of the inner tube being 9,0 x 

0,4 mm, of the outer tube - 13,4 x 0,2 mm). The fuel elements 
were located In alluminiiua spacing lattices. Due to lattice 
pitch variations it was possible to investigate the neutron 
slowing down in the media with different relationship of fuel 
and liquid moderator volumes. 

In this experiment square neutron migration areas values were 
obtained. This values did not differ from the within the 
accuracy of the experiment, the diffusion length values of 
neutrons being negligible small in given multiplying media. 

The value (thermal) is given in Table I 




^56 
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The calculation of square length of slowing down 
by the momentum method in multlgroup approximation 


Neutron tremsfer equation In multigroup approximation with 
elastic and non-elastic slowing down may be giveu as follows /8/: 

i : • i-/ ^ ^ 




where ^ - is group scattering cross section descent from ”y ** 

group; 

r -ia group inelastic scattering cross section from 
ih " group into "y " group; 

[l^—j - ia group sections of neutron slowing down by nuclei 

- source functions; 

- is intergal neutron flux and current within the 

i T 

group; 

ij f/ ” takes into account the neutron multiplication in 
ep.^themal energy region. 

o// - e 

a « . i • e-‘ . g “e 

» ^ ' ly 


\i - i ^-± 




i/^ ue 


Resonance effSets in expression (I) are taken into account 
by adding the following value to leakage sections: 

J'f -- f A-; 

, iy-J i Uj 
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If one adds Fourier transformation to this system (1) and 
then make expansion "PI', we shall obtain a system of relative 
equations for the square length of slewing down determination: 

' / 


J j-i e-j J-i , 

j e=i e.-d 



e>j. 


( 6 ) 



i- 

t 

0 ^/ 


In the expression (3) are space angular moments of the 

neutron distribution function of "y " group in the infinite me- 
dia with the plane isotropic source, 

A program of numerical calculation L\ is made for a com- 
putor/s/.The calculated values of square length of the slowing- 
down in tiydrogeneous media are given in Tables I, II and Pig.1. 

Measurement Results 

" Quasi-homo geneous medi a" . The results of the measurements 
and calculations for pure moderators and quasi-homogeneous 

metal— water media are given in Table I, There is good agreement 
between experimental and calculated data Zij which proves the 
valuability of the calculation method of the neutron slowing do- 
wn as well as miiltigroup constant systems used in calculations 
for the discription of the neutron slowing down process in hyd- 

rogeneous media. 
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Heterogeneous Media , The results of the neutron moderation 
measurements are given In Table II and Pig. I and 2. There is 
great difference In values with the same concentration of 

metal In water. In all cases the measuremented values for 

heterogeneous media exceeded the corresponding values for quasi - 
honogeneous media. The same Is true for diphenyl - metal media 
for which the measured values are much greater than the cal- 

culated ones, obtained by the method applied to homogeneous media. 
This difference can be explained by a sharply expressed anisotropy 
of neutron diffusion in heterogeneous media. 

In the previous papers dealing with fast neutron slowing; down 
In metal -water media was not observed such sharp anisotropy, as 
the investigated media were, actually, quasl-homogeneous /9/. 

The results obtained are evidently of interest for physical 
calculations of reactors, the core of which is sharply heterogene- 
ous. 


II. THERMAL MEUTHON DlPPUSIOg 
Experimental Methods 

Heutron diffusion in hydrocarbons (cyclic and non-cyclic com- 
poimds) has been studied In the medium with a wide range of tempe- 
ratures (from melting to boiling). Por investigation purposes the 
method of pulsed neutron source /II/ has been used. 

Thermal neutron distribution In the medium may be described 
by means of a one-group diffusion equation In which diffusion con- 
stants averaged by the thermal neutrons spectrum are used: 


D(v) i>n v ^ 




Where 


J)(yf)s - diffusion coefficient in the infinite 

medium; 


- neutron absorption rate in the infinite 
medium . 

For a finite medium the changing of thermal neutron density 
in the course of time may be described as a sum of exponentially 
attenuating harmonics. A sufficient period of time having passed, 
all harmonics attenuate, except the main one, and neutron 


distribution is described by the law g 


- cit 
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C/-- r-Vi-M-Jl - (C^-C^) Jl' 


, where decay constant 
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Ji — geometrical buckling; 

— diffusion cooling coeffipient; 

Cy. — non-diffusion correction. 

Description of Experiment 

Measurements were made with the cyclotron of the Academy of 
Science of the USSR.The unit consisted of a cylindrical tank 50 cm 
in diameter, shielded with a B^C layer of 10 mm thickness and with 
a 0,5 iii™ cadimium layer, Into the tank there was inserted a cadmium 
"piston" by means of v/hich the geometrical dimensions of the me- 
dium under investigation could be changed qidLckly and easily. Under 
the tank bottom a block of counters was placed separated and 
shielded from the tank bottom with the help of a cadmium "mask" the 
shape of which can be determined by the function t <?, r) ■ 

This "mask" decreased the influence of high harmonics of the 
neutron density for high geometrical buckling, For heating the me- 
dium an electric furnace was placed on the side surface of the tank. 
For temperature control a thermocouple was placed inside the tank, 
the accuracy of measuring the temperature being +2®C. 

The calculation of the unit and the control tests were carried 
out with water at the temperature of 21®C,For determining the effect 
of mutual arrangement of the cyclotron target and the unit , measure- 
ments were at different distances from and different angles to the 
target. The results remained the same within the experimental error. 
As water shielding and other scalterers were close to the unit, 
but the target was rather far from it, the fast neutron source 
was practically volumetric. 

The liq,\iid was heated to high temperatures and therefore the 
minimum distance between a block of counters and the tank bottom 
could be 2 cm. In connection with this all possible distortions of 
the decay constant were measured because neutrons of various ener- 
gies covered the distance from the tank bottom to the centers in 
different time intervals. 

The results of the experiment showing the effect of the covered 
^is'fcaLnce on the decay constant shows that this effect in negligible 
from JL =0,2A0 cm^.The measurements were for JL which 
were not larger than this quantity. 

The cyclotron had a pulse-repetition rate of 500 cycles and 
the pulse duration of 6-8 jn sec. 

Time analysis of the decay of thermal neutron flux was made 

with the 60-channel time analyzer, the triggering of which 

356 
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was synchronous vrith the cyclotron start-up. The delay of the 
time-analyzer trigij;ering relative to the fast neutron pulse 
varied with-in 100-300 yl< sec. The time channel width was 2-8^ sec. 
The vibration of the base was not more than 0,25 yv sec. 

The method of treatment the results 


where 


.DM, 




The values v;ere calculated 

by the method of least squares. The dependaiace of the diffusion 
coefficient from the temperature of the medixxm for all the inves- 
tigated materials is well discribed by the law: 

(9j 

- is the diffusion coefficient , averaged by the 
spectre of thermal neutrons at the temperature 
of medium T; 

- is the density of the medium at given tempe- 
ratures T,To; 

Q - is the parameter, depending upon the molecule 
structure of the material. 

The results of the treatment of experimental data are given 
in table III. 

It is natural to suppose, that the established neutron spect- 
rum in an infinite medium oveys to the Maxwell distribution: 

where IT • -is the most probable velocity of the 

neutrons, 

Strickly speaking, the spectre of the thermal neutrons obeys 
to the Maxwell distribution only in nonabsorbing media. But 
as the experiment shows one can consider th( spectre as a 
Maxwell one if the absorption is weak. This allows to define the 
value within a small mistake. In this case one uses the 
medim temperature as a parameter. 

Then the equation /^/ may be written as follows: 

4 MTr = ^ K M 

The diffiasion coefficient averaged by the Maxwell spectra 
be written in the form ofs 


can 


DM - 




(ii) 
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Prom the equation (10) and (11) one can get the ratio (11^ 
for defining the differential value of the diffusion coeffi- 
cient of neutrons with the velocity if by the values of the 
diffusion coefficients of neutrons for the most probable veloci- 


ty 






3 

nfa) 


K 


M 

M 


where - is the coefficient of the averaged fvmction, 

depending upon the kind of neutron spectrum 
and of the index (f 

The averaged diffusion coefficient of neutrons by the Maxwell 
spectrum can be defined by the differential value of the diffu- 
sion coefficient of neurons for the most probable velocity. 


DM -DMJ 


Je 




v 


Zro 






- /n(a}J)C\r,) 


Prom the equation one can get the dependence of the 

transport length upon the neutron velocity 

eoid we finally have 

h,(a) Dfy.} 


J>(y)r ■ 


a- / 


( 14 ) 

( 16 ) 


It follows from the abovement ioned, t hat the value of the depen- 
dance of the integral value upon the medi\im tempera- 

ture and the knowledge of the kind of the spectrum of the 
thermal neutrons allows to get the differential value 
from the equation • 

And the differential value (iPo) from the equation 

Prom the equations /16/ and /16^/ it follows 

4 ■ 

^ \fr h)(a-i) s 


(I8j 
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Therefore it is necessary to note in references (e.g.11,13) 
an in accurate ratio of t he diff usion coefficient is used 

W - (^^) 

Got from the equation {^^] has no defin ite phy sical sense. 
This concerns also to the definition of {vj received 

from the experiments with stationary diffusion, In this case 
is an average of the spectrum of thermal neutrons: 

= {2.0] 

■Jf yn(ei-i) 3 \ / 

Then in order to get an information of the transport length 
value, it is necessary to fin d the value of .This can be 

done j.f one defines at different temperatures of the 

medium and the results are treated according the abovermentioned 
method. 

Results 

From the values of the transport length one can calculate 
the transport cross section of the neutron scattering by the 
molecules of the medium, averaged by the neutron spectrum; it is 
also possible to obtain the differential value of the transport 
cross section of moleculs neutrons of the raedim in the case of 
the Maxwell distribution. The latter makes it possible to cal- 
culate the transport cross section Of one hydrogen atom, consi- 
dering the atoms of carbon being free. 

Thus it is possible to determine both the integral value 
of the transport cross section of the neutron scattering by the 
boimd hydrogen and the differential cross section. The latter 
values must correspond to the values of the transport cross sec- 
tion obtained from the differential experiments. 

In order to verify the methods of treatment of the experimen- 
tal results, a comparison of the values of transport cross section 
of neutron scattering by the medium molecules was made, These 
values were obtained both from integral and differential expe- 
riments. The comparison was made for water. The results of the in- 
tegral experiment for water /11/ were treated according to the 
abovementioned method /equation 17/« 

Prom differential experiments on the scattering of neutrons 
with different energies with water COS Q was talcen from /12/ 
and 6y ~ from /G3/ ,and (5^ was calculated, was 

also obtained from /11/ using the expresflinn ("lyl.The results 
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agreed perfectly. Besides, from /11(5)/ s;. were obtained 

using [19) .These values of differ greatly from the re- 

sults of the differential experiment. The difference is of about 
(fig. 5). 

Table III gives the values of the parameter of O and 
by bound hydrogen for the investigated matter at the medium tem- 
perature i8®C , including the data of water and dowtherm /11,1V 
treated accordi.j^ to the abovementioned method, 

H 

Transport cross section of the bound hydrogen D/g, changes 
from 53 bams for d 1 to 100 barns for a 3. There is a 
whole spectrum of intermediate values (fig. 4) depending on the 
structures of compounds. 

Transport cross section of the neutron scattering by the bound 
hydrogen strongly depends on the medium temperature in the case 
of non-cyclic compounds. 
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TAB, III 
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5. Transport cross secti- 
on of the scattering by 
water, received by different 
methods, 1 -differential ex- 
periment; 2-integral experi- 
ment received from the equa- 
tion (17); 3-integral esqperi- 
m«nt received from the 
equation (19). 


slowing dovm of the 
^.I'o+Be; source in metalic + 

water media; 1-"quasihomoee- 

neous" medium Fe+H5®;2-hete- 
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Fig, 4. The dependence of 
upon the parameter a 
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